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1 Introduction 



In the modelling of physical and other phenomena differential equations, be 
they ordinary or partial, scalar or a system, are a common outcome of the 
modelling process. The basic problem becomes the solution of these differen- 
tial equations. One of the fundamental methods of solution relies upon the 
transformation of a given equation (or equations; hereinafter the singular 
will be taken to include the plural where appropriate) to another equation of 
standard form. The transformation may be to an equation of like order or of 
greater or lesser order. In the early days of the solution of differential equa- 
tions at the beginning of the eighteenth century the methods for determining 
suitable transformations were developed very much on an ad hoc basis. With 
the development of symmetry methods, initiated by Lie towards the end of 
the nineteenth century and revived as well as developed ever since, the ad 
hoc methods were replaced by systematic approaches. About the same time 
classes of equations were established as equivalent to certain standard equa- 
tions. In particular the possibility that a given equation could be linearised, 
i.e. transformed to a linear equation, was a most attractive proposition due 
to the special properties of linear differential equations. Already in his thesis 
of 1896 Tresse E3] showed that the most general second order ordinary dif- 



ferential equation which could be transformed to the simplest second order 
equation, videlicet 

X"(T)=0, (1.1) 
by means of the point transformation 

X = F(x,t) T = G(x,t) (1.2) 

has the form 

x + A 3 (x, t)x 3 + A 2 (x, t)x 2 + Ai(sc, t)x + A (x, t) = 0, (1.3) 

where the overdot denotes differentiation with respect to the independent 
variable t whereas the primes refer to T derivatives. This notation is used 
throughout the paper. In terms of the transformation functions F and G the 
functions A, are given by 

A,i.r./ j = /•:,/;,- /••/;,,; u 

A 2 (x, t) = [F xx Gt + 2F xt G x — 2F x G xt — F t G xx ] / J 
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A 1 {x,t) = [F tt G x + 2F xt G t -2F t G xt -F x G tt }/J 

A (x,t) = [FuG t -F t G tt ]/J, (1.4) 



where J(x, t) = F x G t — F t G x 7^ is the Jacobian of the point transformation 
( |1 . 2| ) . As usual subscripts denote partial derivatives. Through the elimina- 
tion of the transformation functions F and G it is found that the coefficient 
functions, A,, % = 0, . . . , 3, must satisfy the conditions 

A lxx - 2A 2xt + 3A 3tt + 6A 3 A 0x + 3A A 3x - 3A 3 A lt - 3AxA 34 
-A 2 A lx + 2A 2 A 2t = (1.5) 

A 2tt - 2A lxt + 3A 0xx - QA A 3t - 3A 3 A 0t + 3A A 2x + 3A 2 A 0x 
+AxA 2t - 2A 1 A lx = 0. (1.6) 

The problem has attracted some interest since the work of Tresse. See & ^(| 
|], Note that the condition corresponding to ( |1.5|) given in [|J omits the 
coefficient 2 in the final term. 

In a practical context one identifies the coefficient functions from the given 
nonlinear equation and substitutes them in the compatibility conditions (|1.5| ) 
and ( |1.6| ) to determine whether the equation is of the correct form. If this be 
the case, the transformation is determined from the solution of the system 
( |1.4| ) and the solution to the nonlinear equation follows immediately. Since 
all second order linear ordinary differential equations are equivalent to (|1.1| ) 
under a point transformation (i.e. of the type ( |1.2| )), these formulae answer 
the question of the class of second order equations linearisable under a point 
transformation. 

The attraction of point transformations is that they preserve Lie point 
symmetries. Since a second order linear equation possesses eight Lie point 
symmetries, which is far in excess of the two required to reduce the equation 
to quadratures, there is no necessity to confine one's interest to point trans- 
formations if the matter of interest is the solution of the equation and not its 
point symmetries. One can then look towards some type of transformation 
which is more general than a point transformation. The advantage of point 
transformations is that they are fairly easy to work with. In looking towards 
a generalisation one wants to keep this property, if possible. 

A convenient form of generalisation is the nonpoint transformation intro- 
duced by Duarte et al [§] which has the form 

X(T) = F{t, x), dT = G(t, x)dt. (1.7) 
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Without a knowledge of the functional form of x(t) the transformation ( |1.7| ) is 
a particular type of nonlocal transformation. This transformation is a gener- 
alisation of a transformation proposed by Sundman [ 21 \ . Since the expression 
'non-point transformation' has a meaning more general than that of (|1.7|) , 
it may be better to refer to this type of transformation as a generalised 
Sundman transformation. 

In their paper Q Duarte, Moreira and Santos derived the most gen- 
eral conditions for which a second order ordinary differential equation may 
be transformed to the free particle equation X" = under the generalised 
Sundman transformation ( TOT ). Euler H] studied the general anharmonic 
oscillator 

x + fx(t)x + f 2 (t)x + h{t)x n = (1.8) 

and derived conditions on the coefficient functions fj for which ( |1.8| ) may be 
linearised under the generalised Sundman transformation ( |1 . 7|) . It follows || 
that ( |1.8| ) may be reduced to the linear equation 

X" + k = 0, ke^\{0} (1.9) 

by the transformation 

X(T) = h{t)x n+ \ dT = ^Lt±f 3 (t)h(tfj x n dt, (1.10) 
where n G Q\{-3, -2, -1, 0, 1} and 

ft(()=/r .)/(» + 3» exp | 2 ^ j' h{p)dp ), (1.11) 

if and only if fx, f 2 and fs satisfy the following condition: 

1 fo n.4-4 f h\ 2 

h 




3)a ,x- (1-12) 

This leads to the invariant (time-dependent first integral) of ( |1.8j ) through 
the first integral of ( |1.9| ), which is 

I(X,X') = X + ±-(Xf. 
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Within the two classes of transformation given by ( |1.2| ) and ( |1.7| ) there 
are subclasses which have some particular interest if one is concerned about 
the type of transformation and preservation of types of symmetry. For ex- 
ample in the case that G x = both classes of transformations reduce to a 



transformation of Kummer-Liouville type [|K], [16[] which has the property of 
preserving symmetries of Cartan form, i.e. fibre-preserving transformations 
[[]]. These transformations are of importance for Hamiltonian Mechanics and 
Quantum Mechanics. 

In the studies of second order equations use was made of the fact that 
every linear second order equation is equivalent under a point transformation 
to (|1.2p . We recall |l5j [p 405] that the number of Lie point symmetries of a 
second order equation can be 0, 1, 2, 3 or 8 and that all linear second order 
equations have eight Lie point symmetries with the Lie algebra sl(2, R). Nat- 
urally any second order equation linearisable under a point transformation 
also has eight Lie point symmetries. In the case of third and higher order 
equations such economy of property does not persist. An n-th order linear 
differential equation can have n + 1, n + 2 or n + 4 Lie point symmetries 
14]] . Consequently, even under point transformations, there are three equiv- 



alence classes of linearisable n-th order equations. In addition there are the 
classes of equations corresponding to other numbers of Lie point symmetries 
or different algebras. 

The manipulations required for the calculations of the classes of equa- 
tions which are equivalent under either of these classes of transformation are 
nontrivial. Although, in principle, the same ideas - indeed more compli- 
cated types of transformation - can be applied to differential equations of all 
orders, the burden of calculation has in the past made the widespread use 
of these methods impracticable. These days with the ready availability of 
symbolic manipulation codes on personal computers of reasonable computa- 
tional power the drudgery of these calculations has been removed or, at least, 
transferred to a third nonvocal party. The time has come to consider the pos- 
sibilities of more complicated transformations of more complicated equations. 
In this paper we report the results for generalised Sundman transformations 
for third order ordinary differential equations. In our calculations we make 
use of the packages Crack written in the computer algebra system Reduce 
and RlF written in Maple. They are described in [p3| and [p 
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We concentrate on the basic third order equation 



X"\T) = 0, (1.13) 

but there is no reason why the ideas presented here cannot be extended to 
the other types of third order equation which constitute the set of equivalence 



classes of third order equations |17|. We give only some examples of more 
general linear third order equations (see Section 3). Note that ( |1.13| ) admits 
the two first integrals 

I X = X\ I 2 =X"X- l -{X'f . (1.14) 

The generalised Sundman transformation ( |1.7|) can then be used to express 
these first integrals to obtain the invariants of the nonlinear equation derived 
by the transformation. 

In the case of scalar third order equations there can be 0, 1, 2, 3, 4, 5, 6 and 
7 Lie point symmetries. When the maximum number of Lie point symmetries 
is found in an equation, this equation has in addition three irreducible contact 
symmetries. The ten symmetries possess the Lie algebra sp(5) jflj. The 
equation ( |1.13[ ) is a representative of this last class. 



We recall that the most general linear differential equation of third order 

is 

x + U(t)x + h{t)x + f 2 {t)x + h(t) = 0. (1.15) 
The condition for ( |1.15| ) to be reduced to 



X"'(T) = (1.16) 



by an invertible point transformation was found by Laguerre in 1898 [[□]] to 

^ l£U + 1 f dU lrf/3 + 2 f3 1 ff + f-0 (117) 

We note in passing that it is possible to transform ( |1.15|) to ( p.,16| ) without 



any conditions on the coefficient functions if contact transformations are used 
either directly or as an equivalent point transformation of the corresponding 
linear first-order system . 

The article is organised as follows: In Section 2 we present the classes of 
equation equivalent to ( |1.13| ) under the generalised Sundman transformation 
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( |1 . 7| ) . In Section 3 we consider a special Sundman transformation and show 
how linearisable ordinary differential equations of second, third and fourth 
order and their invariants can be constructed. Section 4 is devoted to an 
extension of the generalised Sundman transformation, which is related to 
the so-called generalised hodograph transformation introduced recently for 
evolution equations ||. We note that there has been work done recently 
on the problem of transitive fibre-preserving point symmetries of third order 
ordinary differential equations 0, as well as contact transformations and 
local reducibility of an ordinary differential equation to the form ( p.. 13 ) [!]• 
The results presented here are complementary to the results presented in 
those two papers. 

2 Generalised Sundman transformations for 

X'" = 

We turn our attention to the equivalence class of third order nonlinear ordi- 
nary differential equations obtainable from ( |1.16| ), videlicet 



X"\T) = 

by means of the generalised Sundman transformation ( |1 . 7| ) , videlicet 

X(T)=F(t,x), dT = G(t,x)dt, 

where F, G G C 3 and are to be determined for the transformation of (|1.16|) . 
Provided F X G 2 ^ 0, the form of the representative equation of the equivalence 
class is 

x + A 5 (x, t)x + A 4 (x, t)xx + A 3 (x, t)x 3 + A 2 (x, t)x 2 + Ai(x, t)x + A (x, t) = 0, 

(2.1) 

where the functions Aj are related to the transformation functions F and G 
by means of 

■ / \ „ F T t Gf Ft G r 
A4OM) = -4^ + 3- F - 



G F x 



8 



A 3 {x,t) 
A 2 (x,t) 

Ai(x,t) 
A (x,t) 



3 ^ "I - 3—— 



-L X -L X 

G x Gt 
~G~G ~ F n 



G x 
G 
Ft G xx 



F xx G x G 

XX 

~F~~G~G 

n F x t G x 



, F xx Gt 
'~Fx~~G 



G 



xt 



3 

-2 



Ft G xt 

'fZ~g~ 



x G 
F x tt 



Fx 
Gu 
~G 



G 

G x „ F xt G t , „ F t G x G t 

3 ^- 6 ^ +6 f:^g 



(2.2) 



F x G 



Ft Gtt Ft I ( i i . , , , , 

--f -£ - 3-f -4 + ^ + 3-f — ^ 



tft 



Ft (G 



F x G 



F x G 



G 



In order to have the inverse transformation we must establish the compati- 
bility conditions for ( |2.2j ). This is the major thrust of our code. 
The above set of equations has the form 














E x 
E 2 
Ez 
E 4 

E, 

E R 



3F tx G — 3F x Gt — F t G x — F x GAz, 
3F XX G — AlF x G ' x — F X A 4 
F XXX G 3F XX G X G , , -v ... , •, -v ... • 
3F txx G 2 — 6F tx G x G — F t G xx G + 3F t G 2 — 3F xx GtG 
—1F x G tx G + 6F x G t G x — F X G 2 A 2 



F X G XX G + 3F X G 2 — F X G 2 A 3 



-3F ttx G 2 + QFtxGtG + 3F U G X G + 2F t G tx G - 6F t G t G a 
+F X G U G — 3F X G 2 + F X G 2 A 1 
FutG 2 — 3F a GtG — F t GuG + 3F t G 2 — F x G 2 Aq 



(2.3) 
(2.4) 
(2.5) 

(2.6) 

(2.7) 
(2.8) 



In the following simplifications of this system the equation to be replaced in 
each step is multiplied with a nonvanishing factor so that the new system 
after each replacement is still necessary and sufficient. 



E 3 - 


- E 7 : 


= E 2x G — 


G X E 2 — 3E 3 


E 7 - 


+ E 8 : 


= ((5G, - 


GA 4 )E 2 - 3E 7 ) /F x 


E, 


-> E 9 : 


= E 2t G — 


G t E 2 — E4 


Eg ~ 


-> E 10 


:= —G t E 2 


+ E 9 


E\o - 


- E n 


:= (GA 4 - 


- 2G x )Ei + 3E W 
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E 5 - 


-> E12 


:= G^it — Gj-Ei + -E5 


E12 - 


-> E 13 




En ~ 




:= -Ffi? 8 — 



The new system consists of equations Ei, E2, Eq, Eg, Ei 3 , E u and is trans- 
formed to two new functions h and p which are related to F and G through 
the relations 

F(x,t) = p(x,t)h' 1 (x,t) (2.9) 
G(x,t) = h~ 3 / 2 (x,t). (2.10) 

After making all equations free of a denominator and performing three sim- 
plifications 

E13 — > £15 := (Ei 3 — 3h t Ei)/h 
En — > -E16 := (-E14 — 3h x Ei)/h 
E 2 -> £ 17 := (p£ 8 - 3£ 2 )//i 

we introduce new functions A 6 ,A 7 ,A 8 through 

A 6 (x,t) = -6A 5t + 6Ax - 2A 2 5 (2.11) 
A 7 (x,t) = 6A 4t - 6A 2 + 2A4A5 (2.12) 
A 8 (x,t) = -6A 4x + I8A3 - 2A*. (2.13) 

To compactify the display of the resulting system we use the notation 

:= A — Alp^h 

where A is a differential expression in the functions p and h and A\ p ^h is the 
expression after swapping p and h: 






= E 8 - 


= 9h xx - 3h x A A + hA 8 


(2.14) 





= En 


= QPxx - 3p x A 4 + pA 8 


(2.15) 





= Eq - 


= [2h ttt p + 3h tt p t - 2h x pA ] p ^ h 


(2.16) 





= Ei~- 


= [6h tx p - 3h t p x - 2h x pA 5 ] p ^ h 


(2.17) 





= Eiq 


= [18h tx p x - h t pA 8 + h x pA 7 ] p ^ h 


(2.18) 





= Ei 5 


= [Qht x pt - lSh tt p x + 3h t p x A 5 + h x pA 6 ] p ^ h . 


(2.19) 
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Remarkably the system is symmetric under the exchange p «-> h. 



Remark: In view of the symmetry p «-> h and the relations 
we obtain the transformation coefficients 

F{x,t) = F-\x,t) 
G{x,t) = F- 3 / 2 (x,t)G(x,t) 

for the generalised Sundman transformation fjl.TQ . This does not lead to new 
linearisable third order ordinary differential equations, so we do not list here 
any conditions for this transformation. 



jz) and (inn) 



We obtain a first integrability condition by differentiating equation ( p. 17] ) 
with respect to x and substituting h txx , h xx using equation fl2.14|) , substitut- 
ing ptxx,Pxx using equation ( |2.15| ), substituting h tx p x using equation ( [2.1S| ) 
and substituting h tx p using equations (|2.17|) . The result is 

= E 18 := [{h x {\2K u - 12A 5 , - A 7 ) - h t A 8 )p] p ^ h . (2.20) 

Before we give the most general solution we consider two special cases in the 
next two subsections. 



2.1 The case G x = 0, 

i. e. the transformation 

X(T) = F(x,t), dT = G(t)dt, 

where F(x,t) = p(x, t)h~ l {t) and G(t) = h~ 3 / 2 (t). By investigating the case 
G x = (which is equivalent to h x = 0) we cover the case p x = as well 
because of the p <-> h symmetry. 

For h x = we have h t ^ and get E 8 = = A 8 and further E w = = 
A 7 . After substitution of p tx from equation (|2.17|) into equation fl2.19|) the 
resulting system is 

= -E 18 /{12p x ) = A At -A 5x (2.21) 

= E 19 := E 17 /3 = 3p xx - Px A 4 (2.22) 

= E 1 = -Qhp tx - 3h tPx + 2hp x A 5 (2.23) 

= E 20 := (3/i^i - 2hE 15 ) jp x = 36h tt h - 9h 2 t + 2h 2 A 6 (2.24) 

= E 6 = 2h ttt p + 3h tt p t - 3h t p tt - 2hp ttt + 2hp x A . (2.25) 
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For equation ( f2.24| ) to have a solution for h(t) the condition on Aq is 



A 6x = 0. 



(2.26) 



We need to derive one more integrability condition before being able to for- 
mulate a procedure to solve the above system. We reduce the condition 



E, 



6x 







(2.27) 



by substituting ptttx,Pttx,Ptx computed from equation ( |2.23| ), substituting p x 
from equation (|2.22j ) and substituting h tt from equation ( |2.24|) to get 



= (108/i 2 £ 6z - \2QhhuEx + 54h 2 t E l - 3Qhh t E lt + 18h t p x E 20 
-QhhtE^ + 36h 2 E m + 12h 2 E lt A 5 - 9hp x E 20t + 2Ah 2 E x A 
-6hp x E 20 A 5 - 72h i E 19 A + Ah 2 E l Al)/{2h'p x ) 
= 108A te - 36A 5 « - 36A 5t A 5 - 9A 6t + 36A A 4 



-AA\ - 6A 5 A 6 . 



5t 

(2.28) 
(2.29) 



We summarize: The procedure for a given set of expressions Ao, Ai, . . . , A5 
is as follows. 

1. Compute A 6 , A 7 , A 8 from equations (PUD , (pl2|) , (|2~T3|) . 

2. The following set of conditions for Aj is necessary and, as becomes clear 
below, also sufficient for a solution with h x = = G x to exist: 

A 7 = 0, A 8 = 0, A At -A 5x = 0, A 6x = 0, 

108Aox - 36A 5 « - 36A 5f A 5 - 9A 6t + 36A A 4 - 4A 3 - 6A 5 A 6 = 0. 

3. The function h(t) is to be computed from the condition ( |2.24j ) which 
is an ordinary differential equation due to Aq x = and which can be 
written clS db linear equation for /1 3 / 4 : 

24(hV% + h^A, = 0. 
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4. Compute a function q(x,t)(= \og(p x )) from the two equations Q2.22| ) 
and ( p.23|) as a line integral: 



1 A 1 A 

:A 4 , Qt = 77 A 5 



lht 

2 h ' 



«. 3 - - 3 
The existence of q is guaranteed through condition ( |2.21| ). 
5. Compute a function r(x,t) from 

r(x,t) = / exp i)] dx. 



6. The function p(x,t) is computed from 

p(rr, t) = r(x,t) + s(t), 



where the function s(t) is computed from equation ( 2.25|) which after 
the substitution p = r + s takes the form 

2h ttt (r + s) + 3h tt (r + s) t - 3h t (r + s) tt - 2h(r + s) m + 2hr x A = 0. 

This is a linear third order equation for s(t). It is an ordinary differen- 
tial equation with a solution for s(t) because the derivative of the left 
hand side of this equation with respect to x vanishes due to equations 
(IPD , (]pD, (ET23D , (|gg ) and (^29|) . 



7. With the last step we satisfied all conditions (|2.21|) - (|2.25| ) and com- 
puted p and h which gives G and F through equations ( |2.10| ) and 

G{t) = h- 3 / 2 , F{x,t)=ph- 1 . 

2.2 The case F t = 0, 

i.e. the transformation 

X(T) = F(x), dT = G(x, t)dt. 



For all investigations below we can assume h x ^ 0,p x ^ 0. We consider 
ft = (p/h) t = so that 

■ V -. (2.30) 
p 



h f = h- 
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A straightforward substitution of h t from equation ( |2.30| ) into equation ( p. 17] ) 
gives 

9 Pt - 2pA 5 = 0. (2.31) 
Application of this condition on equation ( |2.30| ) yields 

9h t - 2hA 5 = (2.32) 

and both together simplify equations fl2.14j ) - ( ^.19[ ) to 



A = (2.33) 

4A 5x - A 7 = (2.34) 

12A M + 2A 5 2 + 3A 6 = (2.35) 

9p xx -3 Px A 4 + pA s = (2.36) 

Qh xx - 3h x A 4 + hA 8 = 0. (2.37) 

Substitution of h t ,p t into the integrability condition ( |2.20| ) gives 

12A 44 - 12A 5X - A 7 = 

which when simplified with equation ( |2.34|) yields 



3A 4t - AA 5x = 0. (2.38) 
The remaining integrability condition between equations ( |2.31| ) and ( [2.36] ) 



(and equally (|2.32j ) and (|2.37|) ) is computed by the differentiation of Q2.36j ) 



with respect to t and the replacement of Ptxx,Ptx,Pt with ( |2.31| ), p xx with 
(^36|) and A 4t with (|3§. The result is 



6A 5xx - 2A 5x Az + 3A 84 = 0. (2.39) 

We summarize: The procedure for a given set of expressions A , A 1; . . . , A 5 
is as follows. 

1. Compute A 6 , A 7 , A 8 from equations ( gXTJ , (p32| , ( |2TT3|) . 

2. The following set of conditions for A« (i = 0, . . . , 5) is necessary and 
sufficient for a solution with F t = 0, h x ^ to exist: 

A = 0, 4A 5a; -A 7 = 0, 12A 5i + 2A 2 r> + 3A 6 = 0, 

3A 4t - 4A 5x = 0, QAs xx - 2A 52 A 4 + 3A 8t = 0. 
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3. A function u(x, t) is to be computed as 



u(x, t) = exp 



- / A 5 (x,t)dt 



4. The functions p(t,x) and h(t,x) are computed from 

p(x,t) = v(x)u(x,t), h(x,t) = w(x)u(x,t), 

where v(x) and w(x) are to be computed from the two linear second 
order conditions 

9( vu )xx ~ 3{vu) x A 4 + (vu)A 8 = 0, 9(wu) xx - 3(wu) x A 4 + (wu)A 8 = 

which when divided by u are purely ^-dependent and therefore are 
ordinary differential equations for v(x) and w(x) as guaranteed by the 
integrability conditions above. 



5. With the last step we satisfied all conditions (|2.31|) - ( |2.37| ) and com- 
puted p and h which give G and F through equations ( |2.10| ) and (|2.9|): 



G(x,t) = h~ 3/2 , F(x)=ph- 1 . 



2.3 The general conditions for X 



in 





In order to invert the system ( |2.2j ) in general, that is to solve F and G from 
( |2.2|) , we need to consider three different cases. The obvious conditions which 
apply in all cases are p ^ 0, h ^ 0, F X G ^ 0. Recall further that X'" = is 
transformed into (12 . II) , videlicet 



x + A 5 (x, t)x + A A (x, t)xx + A 3 (x, t)x 6 + A 2 (x, t)x + Ai(x, t)x + A (x, t) = 



and A 6 , A 7 , A 8 are defined in (|2.11|) -( 2.13|) ), videlicet 



A 6 (x,t) = —QA 5t + 6A X — 2Ag 
A 7 {x,t) = 6A 4t - 6A 2 + 2A 4 A 5 
A s (x,t) = -6A 4x . + I8A3 - 2A^. 
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In providing the general conditions in this subsection we are not able to 
document each individual step like for the previous two subcases. Instead 
we give only the final result, i.e. we list the conditions on h and p and the 
conditions on the AjS which must be satisfied for a given third order ordinary 
differential equation in order to establish linearisation to X"'(T) = by (|1 . 7|) . 

Below we use the notation 

{A{h) = 0V P := {A(h) = and A(h)\ h ^ p = 0} , 

where A(h) denotes the differential expression in h. Thus {A}h^ p represents 
two differential expression; one in h and the same differential expression but 
with h replaced by p. 

Case I. A 8 ^ 0, -ph x + p x h ^ 0, 2A 8 A 4 - 3A 8x ^ 0: 

The functions h and p for the transformation 

X(T) =p(x,t)h- 1 (x,t), dT = h- 3/2 (x,t)dt 
are to be solved from the following set of linear equations: 



h xx — -hxA^ — -hAf. 
6 y 



{ kt = 9X| ( 108 ^ A4 * As _ 108/ ^ A 5*A 8 - 9h x A 7 A 8 + 2hA 2 8 A 5 - 2hA 8 A 7x 
+16/iA 8 A 5:E A 4 - 4/iA 8 A 8t - 16/iA 8 A 4t A 4 - 2hA 7 A 8x - 24hA 5x A 8x 
+ 24hA 4t A 8x + 2/iA 8 A 7 A 4 > 



/ h*-*p 

The following conditions on A, (i = 1, . . . , 8) are to be satisfied: 

A 4« = 7^T2 (-24A 8x A 2 5x - 24A 8x A% - 2A 8:c A 7 A5 :c + 48A &c A 4t A 5:c 

+2A 8:[ .A 7 A 44 + 16A 8 A^A 4 - 2A 8 A 4t A 7 A 4 + 16A 8 A§ X A 4 - 2A 8 A 5ir A 7:c 
-AA 2 8 A 5 A 4t - Af A&, + AA 2 8 A 5 A 5x + 2A 8 A 4t A 7x + 2A 8 A 5x .A 7 A 4 
-16A 8 A 8t A 5x . - 32A 8 A 5x A 4t A 4 - A 8 A 8t A 7 + 12AgA 5x . 4 + 16A 8 A 8i A 4t ) 
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A 5 « = ^| (-144A 8a; A5 ;E A 7 A4 t + 38A 8 A^A 4 A 5;E - 96A 8 A 7 A 4 A^ 

+6A 8 A 7;E A 4t A 7 - 96A 8 A^A 4 A 7 - 288A 8 A 7;E A 4 i A 5x - GAgA^A^ 
-84A2A 5 A 7 A 4t - 3456A 8 A^A 4 A 5;C + 60A 8 A 7 A 5:c A 8t - 1440A 8 A 4t A 5:c A 8t 
-38A 8 A?A 4 A 4 * + 84A2A 7 A 5 A 5x - 60A 8 A 7 A 4t A 8t + 192A 8 A 7 A 4 A 4t A 5x 
+3456A 8 A^A 4 A 4t - 54A 8 xA?A 5x + 5184A te A 5a .A^ + 72A 8x A| c A 7 
+72A 8;C A^A 7 + 144A 8 A* t A 7ie + 1152A 8 A^A 4 + 12A*A?A 5 + ^A^A^ 
-252A2A 4t A 7t - 216A 4t A2A 6x + 33A 7 A|A 6;C - 1152A 8 A^A 4 - 12A|A 6 A 5 
+72A A3A 4 - 72A^A 5 A 5t - 5184A 8;C A 44 A^ - 5A 8 A?A W + 720A 8 A^A 8t 
+36A^A 7t A 7 + 4A 8 A3A 4 + 54A 8;c A 4t A 7 ' + 216A 5;c A2A 6:c + 720A 8 A^A 8t 
-AA s A 7x A 2 7 + 252A2A 7t A 5:c - 4A 8:c A 7 - 1728A 8x A^ + 21QA 3 s A 0x 
-18A*A 6t + 1728A 8x AL - 8A|A*) 

A 7tt = — ^ (2A'A 7 A 5 A 7x + 252A2A 5;r A 7 A 4t + 3A?A 8t A 8 A 4 - 3A 8 A 7a; A 7 A 8t 

-3A^A 6:c A 7 A 4 - 36A 8:c A 8 A 6a A 4t + 48A 2 A 4t A 4 A 7t + 10A^A 8t A 5 A 7 
-6A^A 7t A 4 A 7 + SeAfeAgA^As^ + 36A 8:c A 7 A 4t A 8t + 6A 8:c A 8 A 7t A 7 
+3A 8:c A 8 A 6;c A 7 - 24A2A 6:c A 5;c A 4 + 24A2A 6x A 4t A 4 + 2A 8x A 8 A2A 5 
+72A 8x A 8 A 7t A 5x - 72A 8:c A 8 A 4t A 7t - 36A 8:c A 7 A5 ;c A 8t - 48A 8 ! A 7t A 4 A5 :c 
-2A2A^A 4 A 5 + 24A 7 A 5:c A 8t A 8 A 4 + 16A2A 5 A 7 A 4t A 4 - 24A 8x A 8 A 5 A 7 A 4t 
-24A 7 A 4t A 8t A 8 A 4 + 24A 8:c A 8 A 7 A 5 A 5:c - 16A|A 7 A 5 A 5:c A 4 - 18A|A 6:ct 
-126A^\4A 7 - 6A*A 7 A 5i + 2^^ - 126A^A 7 - 648A2a 5:c A^ 
+648A2A 4t A^ + 6A^A 4t + 6A2A 74 A 7:c + 30A2A 7t A 8t - 6A 3 8 A 5x Al 
-9A*A 6 A 5:c + 9A^A 6 A 4t - 6A 8 A 7 A^ - 21A|A?A 5a: - 3A 8x A 2 7 A 8t 
-36A 8 A 5;r A5 i - 4A^A 7 A2 - 12A!A 6x A 5 - 18A3A 7t A 5 + 24A 8 ! A 6:c A 8t 
+3A2A 6:c A 7;c + 36A^A 5t A 4t - 18A A^ + 216A*A^ - A*A? - 216A 2 AL) 

A 8U = ^3 (9A2A 7:ra A 7 - 108A2A 7ra A 4t + 1296A 8 A| B A 8sa! + 9A 8 A 2 7 A 8xx 

-16A 2 A 8t A 7 A 4 + 12A 8 A 8t A 8a; A 7 + 108A 2 A 7 ^ A 5* + 24A 8 A 4t A 4 A 7 A 8x 
-24A 8 A 5x A 4 A 7 A 8:c - 288A 8 A 5x A 4 A 4t A 8x - 192A|A 5a .A2A 4 i + 8A3A 5 A 7 A 4 
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-144A 5:c A 8;c A 8 A 7;c + l72Sklk At A Ax X bx + \2klX lx k u X A + 144A 8 A^A 4 A 8s 
-5A2A 8:c A 8 A4 - 12A2A 7:c A 5:c A4 - 5A2A 7:c A4A 7 + 144A 8 A^A 4 A 8:c 
+144A 4t A 8x A 8 A 7x + 180A2A 4t A 4:c A 7 - 12A 2 A 5 A 7 A 8;C - 6A 7 A 8:c A 8 A 7:c 
+40A2A 7 A2A 5x - 180A2A 7 A 4x A 5:c - 40A 2 A 7 A 4 : A 44 - 2592A 8 A 5:c A 8;c;c A 44 
-216A 8 A 7 A 8;c;c A 4t + 216A 8 A 7 A 8:c:c A 5:c + 1296A 8 A^A 8!M; + l%k\h Qxx 
+2AA\A 7t A A - 36A2A 7t A 8:c - 8QAA 2 8 Al x A Ax + 2AA 2 8 A st A 7x - 216A 7 A^A 5x 
+lQ56A 3 8 A 5x A 4t + 159A3A 4t A 7 + 96X^^1 + 4A 2 8 A 2 A 2 4 - 144A3A 5:c A 7 
+3456A 4t A^A 5:c + 216A 7 A^A 4i - 12A^A 8t A 5 + 18A3A 6:c A 4 - 3QA 2 8 A 6x A 8x 
+9QA 2 s Al t Al - 8QAA 2 8 A 2 4t A Ax - 9A 2 8 A 2 A 4x + 6A*A 6 - 1728A^A^ - 6A 2 7 A^ 
+4A^A 2 + 24A*A 5i - QA 3 8 A 2 - 792AiA|,. + 48A^A^ -864A|A^ - 1728A^A 

A Axt = t^t- (2A 8 A 7x + 2A 8 A 5x A 4 + 18A 8 A 5a;:!; + A 8 A 8t - 2A 8 A 4i A 4 
-A 7 A 8x - 12A 5;E A 8;E + 12A 4t A 8;c ) 

Ajxt = 3W ( 4A s A 8* A 7 A 4 - 12A 8 A 8t A 8a; A 7 - 72A 8 A 4t A 4 A 7 A 8a; - 4A^A 5 A 7 A 4 

+72A 8 A 5a; A 4 A 7 A 8;E + 72A 5a; A 8:c A 8 A 7;E + 48A 8 ; A 7;c A 4t A 4 + 8A?A 8a; A 8 A 4 
-48A|A 7;c A5 :c A 4 - 4A|A 7;C A 4 A 7 - 72A 4 tA 8:c A 8 A 7:c + 12AlA 5 A 7 A8s 
-16A2A 7 A2A 5x + 16A 2 A 7 A2A 4t - 36AlA 6xx - 12A|A 7t A4 + 36A2A 7t A to 
+12A|A 8t A 7x - 72A 7 A^A 5x - 216A|A 5:c A 4t - 66A|A 4t A 7 - 2A 8 ! A;;A 4 ! 
+54A 8 i A 5x A 7 + 72A 7 A| c A 4t - 12Af A^A* + 36A 2 A 6;C A 8:!; - 12A^A 5 A 7:c 
+QA 2 8 A 2 X - 3A^A 6 - QA 2 A 2 8x - 2A$A| - 12A|A« + 3A 3 A 2 
+72AIAL + 144A^y 

Asxt = ~ IsXf (" 8A 8 A 4* A 4 - 12A 7 AL - 144A5XAL + 144A 4t A| c 
+9A 2 8 A 7xx - 36A|A 5;c - A^Ai - 9A|A 7 A 4:C + 12A 8 A 8:c A 5;c A 4 
-12A 8 A 8:c A 4t A 4 + 5A 8 A 8x A 7 A 4 + SA^A 2 + AA 2 8 A st A 4 + 9A 8 A 7 A 8ra 
-24A 8 A 8t A 8;c - 108A 8 A 4t A 8ra + 108A 8 A 5:c A 8;ca; + 27A^A 4t 
-12A 8x A 8 A 7:c + 72A 8 ; A 4t A 4:c - 72A2A 5;c A 4:c ) . 
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Case II. A 8 ^ 0, -ph x + hp x ^ 0, 2A 8 A 4 - 3A 8x = 0: 

The functions h and p for the transformation 

X(T) =p(x,t)h~ 1 (x,t), dT = h- 3/2 (x,t)dt 
are to be solved from the following set of linear equations: 

{hxx = -h x A A - ^A 8 | 

\h t = -1— (324h x A 4t - 324h x A 5x - 27h x A 7 + 6/iA 8 A 5 
I 27A 8 



-QhA 7x - 12hA 8t + 2/iA 7 A 4 )l 



h^p 

The following conditions on A, (i — 1, . . . , 8) are to be satisfied: 



13 4 22 2 10 

■— A 7 A 8t - — A 5x A 7 A 4 - —A 7x A At + -A^A 4 + 



2 12 

+2A At A^ xx — — A 7 A 7x — 2A 5x A 5xx + -A 4;r A 6x + A 4 A 6ra — -A 7 A 5xx 

4 2 82 79 22 

--A 8 A 6x - -A 2 4 A 6x + —A 8t A u - —A 8t A 5x + — A 7 A 4 A 4t 

2 2 
+— A 2 7 A 4 - -A 4t A 5x A 4 

A 4 « = — ^- (-48A 8t A 5a; + 12A 8 A 5 A 5;E - 6A 5;E A 7;C + 6A 7;E A 4t - 3A 8 A 6a; 
+36A 8 A 5a;t - 2A 7 A 4 A 4t - 3A 7 A 8t - 12A 8 A 5 A 4t + 48A 8t A 4t + 2A 5;E A 7 A 4 ) 



1 

216A§ 

2 A A 1 A A \ A A2 1 A A A 2 



A s« = 777^71 (288A 7 A 4 A 4t A 5a; + 252A 8 A 7 A 5a; A5 + 99A 8 A 6:c A 7 



-648A 8 A 6:c A 4i - 252A 8 A 5 A 7 A 4t + 6A^A 4 A 5:c - 144A 7 A 4 A^ - 144A| t A 4 A 7 
+648A 8 A 6:c A5 ;c - 36A2a 6 A 5 - 864A 7x A 4t A 5:c - lSA^A^ + 36A 8 A?A 5 
+216A A2A 4 - 6A2A 4 A 4t - 216A2A 5 A 5i + 180A 8t A 7 A 5;c - 180A 8t A 7 A 4t 
-756A 8 A 7t A 4t + 18A 7;c A 4t A 7 + 108A 8 A 7 A 7t - 4320A 8t A 4t A 5;c - 54A|A 6 i 
+756A 8 A 7t A 5:c + 2160A 8 *A^ + 648A2A * - 12A 7x A? + 4A^A 4 + 432A| t A 7a 



19 



+432A 7a; A| c - 24AjAl + 2160A 8t A| E - 15A 8t A^) 

A 7tt = -g^a (6A 8 A 7t A 4 A 7 + 2A 8 A?A 4 A 5 - 30A 8 A 5 A 7 A 8t + 5AA 2 8 A 6xt 

-756A 8 A 4t A 5:c A 7 - 6A 8 A 7:c A5A 7 + 3A 8 A 6;C A 7 A 4 + 108A2A 5:c A 5t 
+Q3A s A 2 7 A 5x + 18A^A 7 A 5t - 18A 8 A 7:c A 7t + 18A"A 5:c A^ - 63A 8 A 4t A? 
+378A 8 A^A 7 + 1944A 8 A 5:r A' t + 9A 7x A 7 A 8t - 3A?A 8t A 4 + 54A2A 7t A 5 
+378A 8 A^A 7 - 108A2A 5t A 4i - 90A 8 A 7t A 8t + 12A^A 7 A2 - 1944A 8 A 4t A^ 
-72A 8 A 6:c A 8 i + 36A^A 6;c A 5 - 9A 8 A 6:c A 7:c + 27A 5;c A2A 6 - 27A 4t A^A 6 
-18Aj^A 4i + 648A 8 A^ + 54A A^ + 3A 8 A^ - 648A 8 A^ + 18A 7 A^ t ) 

A 8U = t^t- (24A^ t - 6A 8 A 8t A 5 + 60A 8 A 4t A 7 - 54A 8 A 7 A 5a; + 378A 8 A 4t A 5;E 

loA 8 v 

-4A 8t A 7 A 4 - 3A 8 A? + 12A"A 5t + 2A 2 8 A 2 5 + 12A 8t A 7:c + 3A^A 6 - 198A 8 A|, 
-180A 8 A^ + 9A 8 A 6ra - 3A 4 A 8 A 6:C ) 

A 4a;t = g^rr — g^-5a;A 4 + A 5xx + Yg^ 8t + 7j A 4 tA 4 — — A 7 A 4 

A 7xt = - (12A 7:c A 4 A 7 + 36A 8 A 5 A 7a; - 162A 8 A 7 A 5:c + 648A 8 A 4t A 5x 

-12A 8 A 5 A 7 A 4 - 2A?A^ + 9A^A 6 + 6A 2 8 A 2 5 - 18A^ + 36A^A 5i - 9A 8 A? 
+12A 8t A 7 A 4 - 36A 4 A 8 A 6:c + 108A 8 A to - 216A 8 A^ - 36A 8t A 7:c 
-432A 8 A^ - 36A 8 A 7t A 4 + 198A 8 A 4 t A 7 ) 

2 13 1 

A 7xx = -- A 7 A 2 4 - —A 8 A 4t + AA s A 5x + A 7x A 4 + -A 7 A Ax . 

Case III. A 8 = 0, A 7 ^ 0, -ph x + hp x ^ 0, p x ^ 0: 

The functions h and p in the transformation 

X(T) =p(x,t)h- 1 (x,t), dT = h- 3/2 (x,t)dt 



20 



are to be solved from the following set of linear equations: 



h xx — -^h x A^ 



h^p 



2 

(fx 



{ 

{ kt = 9Xf (" 12/lA 6- A 7 - ^hA 7t A 2 - 2/iA 5 A 3 - 129Qh x A 7t A 6x - 432/^A 
+72h x A 2 A 5t - 8Q4h x A 2 7t - lUh x A 7 A 6x A 5 - U4h x A 7 A 5 A 7t + 432h x A 7 A 6xt 
+A32h x A 7 A 7tt - 12h x A 2 A 2 5 - 18h x A 2 A 6 ) 1 . 

The following conditions on Aj (i = 1, . . . , 8) are to be satisfied: 



A 5 «t = ^4 (-yA?A M A 7 « + — A^A'A^ - yA 3 A 6;ct A 5t - —A 2 A m A 

+^A 7 A 6xt A| B + -^A 7 A? t A 7ti + ^AyA 6 A 7ti + ^^A z 7 A &xt A\ 

rj 166 344 6 

+-A 3 A 6;rf A 6 + — A 7 A 7tt A^ + — A 7 A 2 7t A 6xt - -A*A 5tt A 5 
y to to o 

+^-A 3 7 A 5tt A 7t + ^AlA &x A btt - yA?A? tt - ^A 7 A 6tt - yA?A^ t 
+3A 7 A ^ - 64A 7 ,A 6 , - 1a 7 A - ^A 7 A 2 6 _ ^ A 2 t A 2 6x - ^A 7 A^ 

"P^Ai - ^ALA W + |A?A 5t AL - ^A 7 A^A 6 + ^A?A£A« 

44 6 12 194 

+ — A 7 A 7 t A 6x A 5t - -A z 7 AqA^x - — A 7 A A 4 A 7t + — A 7 A 5t A 6:c A5 
lo o o loo 

86. t2 . * 67 t2t . . 368 . 4 56 . 4 . 4 . . 
— 9~ ~~ 30 A 7 A 6- A ^A 6 - "Y^A 7i - — A^ + A 7 A iA 4 

-^-A z 7 A^ x A 7t + ^A^A ^A 5 - ^A^AzxAsx - ^A 2 A\A\ X + A^A A 5a; 

1 97 7 353 

+-A^A A 4 A 5 + — — A 7 A 6 A 7t A 5 - — A 3 A 5 A 6;c A 6 + — A 7 A 5 A 5t A 7i 

O 54U A(\3 loo 

' 18 -A 7 A^A 6;E A5 - ^A^A^A;? - ^A 7 A^A 5 - ^A?A^A 6 



6xt 



9 11 27 1 9 11 45 

+^A 3 7 A 6t A 7t - ^A^A 5t A2 - Y^A 7 A 5t A 6 - ^A 4 7 A 6t A 5 - yA 7 A^A 5 

-j§A 7 ALA 6 + ^A 6t A 6x + ^A?AijA w - ^A 2 A 2 7t A 2 5 + -^A 3 7 A 3 5 A 6x 
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118 A 2 A A A 148 A2 A A A 148 »AA A 

+— A?A 7t A 6:rt A 5 + — A?A 7t A 5 A 7tt + _A 7 A 7t A 6;r A 6;ri 

1 78 AAA A 124 . 2 A A A 154 a 2 A A A 

+— A 7 A 7t A ra A 6:c + — A 7 A 6:rt A 6:r A 5 + — A 7 A 7tt A 6:r A 5 

A-7ttt = To ( —^T^&x^&xt — ~^ 7 ^&xt^b —MxMtMtt — A 7 A 7tt A 6x 

A 7 V 5 5 5 5 

28 14 3 

-—A 7 A 7t A 6xt - -A 2 A 5 A 7tt + — A 7 A 5tt + A 2 A extt + 6A 7t A^ - — A 3 A 0x 

-^A 7 A 6t + ^A 2 7t A 6x + -A 3 6x + —A 3 7t + ^A?A 6 A 5 + ^A 7 A 5t A 5 

1 1 3 2 11 

-— A 7 A A 4 - — A 2 A 2 5 A 6x + -A 7 A 7t A 5 + —A 7 A 5 A 2 6x + — A 7 A 7t A 6a; A5 

11 7 13 

+— A 7 Aj? + —A 2 A 6x A 6 - —A 2 A 6x A 5t + — A 7 A 6 A 7t - -A 7 A 7t A 5t 

A 5xtt = ^ (~ 3 AnUx + ^A?A 6 - ^A^ - ^A?A 5t - h 7 A ext 

19 8 4 2 7 

-— A 7 A 7 « + -A 2 7t + -A 2 X - — A 7 A 6x A 5 + — A 7 A 5 A 7t + A 7 A A 4x 

-A 7 A 5t A 5x + 3A 7 A 0xx - ^A 7 A 6 A 5x - A 7 A 5 A 5xt - ^A 7 A 2 5 A 5x + A 7 A 4 A ^ 
Aq xx = —-A 7 A 5x — ^A 7 + -A 4 A 6a , 
A 4t = A 5x + ^A 7 
A 7x = ^A 7 A 4 . 

Remark: From our general calculations two more cases were derived which 
are not listed in this subsection, namely the case A 7 = 0, A 8 = 0, A 6x = 0, 
—ph x + hp x 7^ 0, p x 7^ and the case A 7 = 0, A 8 = 0, A 6x = ; — ph x + hp x ^ 
0, p x = 0. Both these cases provide identical conditions on the AjS already 
given in subsection 2.1., i.e. for G x = 0. This means that for these ordinary 
differential equations (specified by these Ais) we already have a simpler lin- 
earisation procedure with p x = 0. Recall that our aim is not to provide all 
linearisation procedures for a given ordinary differential equation but rather 
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to find all linearisable equations allowed by 



We end this subsection with an example of a third order ordinary dif- 
ferential equation which we construct to be linearisable by the generalised 
Sundman transformation. We do not claim any physical relevance for this 
equation. 

Example 2.1: The equation 

3xe x — e x .. 3xte x — Ate x — 4 . .. x 2 te x — 3xte x + 3te x + 3 . 3 

r*> I rf I nn rr> I T* 

x(l + te x ) x(l+te x ) x 2 (l + te x ) 



3x 2 e x — 6xe x + 3e x 
x 2 (l + te x 



^ ^1—^ * = ( 2 - 4 °) 



may be reduced to X'" = by the transformation 
X{T) = te x + x, dT = xdt, 

that is 

h(x, t) = x~ 2/3 , p(x, t) = tx~ 2/3 e x + 



x 



1/3 



By ( |1.14|) and the transformation given here two time dependent first inte- 



grals follow for Q2.40| ). This example corresponds to Case I above. 



2.4 On the computations 

Some remarks about the computations with computer algebra are in order. 
The computational part of this section consists in determining all solutions 
F(x,t) and G(x,t) of the system (|2.2|) , i.e. for each class of solutions a set 
of conditions for A ,...,A 5 and a list of instructions of how to compute 
F and G. Both are obtained by investigating integrability conditions, like 
d x F tX x = d t F xxx with F txx and F xxx being replaced by corresponding ex- 
pressions from two equations. An algorithm for completing this task in a 
finite and systematic way is the well-known (Pseudo-) Differential Grobner 
Basis algorithm. To apply it one has to specify a total ordering of all par- 
tial derivatives of all functions by firstly specifying a lexicographical ordering 
between all functions and also between all variables. In addition one can 
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require that the total derivative of a function should have a higher priority 
than lexicographical ordering which was the case in these computations. 

One of the programs that has this algorithm implemented is the package 
RlF by Allan Wittkopf and Greg Reid which was used for the general cases 
I - III. Its strength is that it is very efficient, especially for larger nonlinear 
problems. It can be downloaded from the World Wide Web site with URL 
http : //www. cecm. sf u. ca/~wittkopf /rif .html. 

For the two special cases G x = and F t = the package Crack of TW 
was used. With its possibility to record the 'history' of any derived equation, 
i.e. how it results from other equations, and with its optional fully interactive 
mode of operation it provided a convenient way for simplifying system (|2.2j ) to 
system ( |2.14| ) - ( |2.19| ). This made it easily possible to prove all conclusions 
for these two cases as one could explicitly describe how each condition is 
generated. For the harder more general cases computed by RlF tracing the 
history of equations would be too costly, so one has to trust the computation, 
although the remaining equations for F and G and the conditions for the A, 
have been tested with several third order linearisable ordinary differential 
equations available to the authors. The program Crack can be downloaded 
from http : / /lie .math . brocku . ca/twolf /home/ crack . html. 



3 A special Sundman transformation 

In this section we consider a Sundman transformation of a special form and 
construct several examples of linearisable ordinary differential equations of 
second, third and fourth order. Some well-known equations are shown to be 
included in this construction. Note that we allow here linearisation to a more 
general linear form than X"'(T) = 0, whereby the transformation functions 



F and G in the generalised Sundman transformation (177) are special 



We consider the following special form of the Sundman transformation: 

x{T{t)) = x {ty 

(3-1) 

dT(t) = x(t) n dt, 
where p, n G Q\{0}. The first three prolongations are 
X' = px p - n ' x x 
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X" = p(p - n - l)x p ~ 2n ~ 2 x 2 + px p - r 



? p-2n-l; 

X'" = p{p - n - l)(p - 2n - 2)x p - 3n - 3 x 3 + p(3p -An- 3)x p -' 3n - 2 xx 



+px p - 3n - 1 x . 
We consider several examples. 

Example 3.1: Let 

X" = J , (3.2) 



where Jo is an arbitrary constant. By the Sundman transformation (|3.1|) we 
obtain in terms of the coordinates (x, t) the following equation: 

p{p - n - l)x p - 2n - 2 x 2 + px p - 2n - x x = J . (3.3) 

Equation ( |3.3| ) can thus be linearised by ( |3.1| ). The first integral for ( |3~2] ) is 

I 1 (X,X')=X-^-(X') 2 

so that a first integral for ( |3.3|) takes the form 

2 

h{x,x) = x p -^-x 2p - 2n - 2 x 2 . 
2 Jo 

As a special case we consider p = n + 1. Then ( |3.3|) simplifies to 

x = -^x n . (3.4) 
n + 1 K J 



Solving (fO) for J and considering this Jo this as a first integral of the third 
order equation which results upon differentiation J with respect to t, we 
obtain 

x — nx~ l xx = 0. (3.5) 
Thus ( p\5[ ) can be linearised by 

X = x n+1 , dT = x n dt (3.6) 
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to the equation 

X"'{T) = 0. 



Note that ( |3.5| ) is a special case of the generalised Sundman transformation 
(|1.7|) with F t = 0, derived in general in subsection 2.2. Equation ( |3.5|) admits 
two first integrals, namely Jq (solved from (|3.4|)) as well as 



h(X,X',X") = XX" - i(X') 2 



expressed in the coordinates (x, t) by ( |3.6| ). 
Example 3.2: Let 

X"'(T) = I , (3.7) 

where Jo is an arbitrary constant. By the Sundman transformation (|3.1|) we 
obtain in terms of the coordinates (x, t) the following equation: 

p(p-n- l)(p - 2n - 2)x p ~ 3n ~ 3 x 3 + p(3p -An- 3)x p ~ 3n ~ 2 xx 

+px P-3n-l = j o ( 3 _g) 

Equation ( PTS| ) can thus be linearised by (|3TT|). Two first integrals for ( |3.7j ) 
are 

i 1 (x',x") = x'-±-(xy 

I 2 {X, X', X") =X + \l 2 (X"f - I^X'X" 

so that two first integrals for (|3.8|) are obtained by expressing the above first 
integrals in (x, t) variables using the corresponding Sundman transformation 

As a special case we let p — — 1 and n = —3/2. Equation ( p.8|) then 
simplifies to 

x = -I x- Fj/2 (3.9) 
and it follows that (|3.9| ) can be linearised by 

X = x~\ dT = x~ 3/2 dt (3.10) 
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to the equation 

X"'{T) = I . 



Solving ( |3.9| ) for I and considering I as a first integral of the fourth order 
equation which results upon differentiation J to t, we obtain 

2xx {4) + hx x = 0. (3.11) 

Equation ( |3.11|) plays an important role in the symmetry classification of the 
anharmonic oscillator (12| . By the above construction ( |3.11| ) can be linearised 
by (|3.10|) to the equation 

= 0. 

Equation (|3.11|) admits three first integrals, namely Iq (solved from ( p.9|) ), 
and 

h(X',X",X'") = X'X'" - i(X") 2 

I 2 (X, X', X", X'") = X {X"'f + i(X") 3 - X'X"X"' 



which are to be expressed in the coordinates (x,t) by ( |3.10| ). 
Example 3.3: We consider 

x — x~ l xx — Aax 2 x = 0, (3-12) 



where a is an arbitrary constant. Equation ( |3.12| ) is of Rikitake type JD3 
By the Sundman transformation (|3.1| ) with p = 2 and n — 1, i.e. 



X = x 2 , dT = xdt, (3.13) 

it follows that (|3.12j ) linearises to 

X m -4aX' = 0. (3.14) 

More examples of linearisable equations using the special Sundman trans- 
formation (|3.1|) are given in pfl. 
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4 An extended Sundman transformation 



In this final section we propose a further extension of the generalised Sund- 
man transformation ( |1.7|) . 

Consider the following exact one-form: 

dT(t, x(t)) = G\(t, x)dt + G2(t, x, x, x, . . .)dx, 

where G\ is a smooth function of x and t, whereas G2 is a smooth function 
of x, t and a finite number of derivatives of x with respect to t. Here d is 
the exterior derivative and by the Lemma of Poincare it is known that 

d{dT) ee 0. 

To find the relations between Gi and G2 we calculate 

d(dT(x, t)) = dGi A dt + dG2 A dx 

dG\ <9(jr2 9G2 . <9(jr2 7.. 

= — — dx A dt H — - — dt A dx H — - — ax A ax H — - — ax A ax + ■ ■ • 

ax dt ox ox, 

dG\ OG2 <9Cx2 .. dG2 ... 

= — — ax A dt H — — at A ax H — - — xdt A ax H — - — x at A ax + • • • 

ax at ax ax 

EE 0. (4.1) 



It follows that 



Also 



dGx _ dG 2 ..dG 2 ...dG 2 

dx dt dx d'x 



dT 

— = G 1 (t,x) (4.2) 
dT 

G 2 (t,x, x, x, . . .). (4.3) 



dx 

We now propose the transformation 



X(T(t,x)) = F(x,t) 

dT(t, x) = Gi(t, x)dt + G2(t, x, x, x, . . .)dx, (4.4) 
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with 



dG 1 dG 2 dG 2 ...dG 2 

+ x——. — h x — — + 



(4.5) 



dx dt dx d'x 

and name this the extended Sundman transformation. Note that with 
G 2 = 0, ( fOj) takes the form of the generalised Sundman transformation 

c 



Next we derive the first two prolongations of ( |4.4j ): 

X'(G 1 + xG 2 ) =F t + xF x 
X"(G 1 + xG 2 f 



+X' Git + iG ix + x [G 2t + xG 2x + xG 2x + x G 2x + 
= F tt + 2xF xt + x 2 F xx + xF x . 
Substituting G 2t from ( |4.5| ) in ( ^4.7| ) we obtain 



X"(G l + xG 2 f + X' G lt + 2xG lx + x 2 G 2x + xG 2 



-\- 2,xF x t ~\~ x F xx -\- xF s 



(4.6) 

xG 2 
(4.7) 



(4.8) 



We give two examples. 

Example 5.1: Consider the first order ordinary differential equation 

X' = X 

and let 

F(t) = t, Gi(x,t)=xt 
so that x(t) = t~ 2 . By ( flBD we obtain 



(4.9) 



Insert G 2 given by ( [4.9|) and G\ = xt in ( f4.5| ) to obtain the second order 
ordinary differential equation 



1 \1 x f 1 

t 2 J X X 2 \t 
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which may be written in the form 



x I xt] + x 2 t + x ( — + x ] = 0. 



t J \t 2 

Obviously x = t~ 2 satisfies this second order ordinary differential equation. 
The transformation is 

x = (XT 1 e- T , t = e T . 

Example 5.2: Consider the second order ordinary differential equation 

X"(T) = 

and let 

F(t) = t, Gi(x,t)=xt 
so that x(t) = t~ l . By ( |4.8| ) we obtain 

x + 2xt + xG 2 + x 2 G 2x = 0, (4.10) 



where G 2 = G 2 {t ). We now solve G 2 from ( |4.10|) by integration. For 

example, a solution is 

2 

G 2 (t,x,x,x) = -* + (*) -2?. (4.11) 

x \xj X 



Insertion of G 2 , given by ( 4.11 ), and G\ = xt into ( |4.5|) gives the third order 



ordinary differential equation 

x'x x + Itxx x — 2x 2 x — 3tx 3 = 
for which x = t _1 is a solution. The transformation is 

x = (X'y l T~\ t = T. 

We note that the linearisations which follow from the extended Sundman 
transformation fl4.4p given in the examples above are 'downwards', i.e. to a 
lower order ordinary differential equation. 
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